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Abstract 

Hadron collider cross sections for ti production and di-lepton, single-lepton and all-jet decays with 
up to 2 additional jets are calculated using complete LO matrix elements with 6-, 7- and 8-particle 
final states. The fixed- width, complex-mass and overall- factor schemes (FWS, CMS & OFS) are 
employed and the quality of narrow-width and double-pole approximations (NWA & DPA) is 
investigated for inclusive production and suppressed backgrounds to new particle searches. NWA 
and DPA cross sections differ by 1% or less. The inclusion of sub- and non-resonant amplitudes 
effects a cross section increase of 5-8% at pp supercolliders, but only minor changes at the Tevatron. 
On-shell tt/Wtb backgrounds for the H — > WW decay in weak boson fusion, the hadronic r decay 
of a heavy and the <fi — > hh — > rrbb radion decay at the LHC are updated, with corrections 
ranging from 3% to 30%. FWS and CMS cross sections are uniformly consistent, but OFS cross 
sections are up to 6% smaller for some backgrounds. 
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I. INTRODUCTION 



Since the discovery of the top quark in 1995 0, the capabilities of ongoing and forth- 
coming collider experiments have improved significantly. Consequently, ti production will 
be abundant and studied intensely as a signal at Fermilab's Tevatron collider and even more 
so at CERN's Large Hadron Collider (LHC). With a decay width T t of about 1.5 GeV the 
top quark decays too rapidly to be observed directly, and is instead identified through char- 
acteristic detector signatures with isolated leptons and jets. These signatures would also be 
observed in the production of various hypothetical particles, so that top production consti- 
tutes an important background for many new particle searches. In light of the changed role 
that top production will play in the near future, the quality of corresponding theoretical 
predictions needs to be reviewed. 

As is well known, a general, systematic and "natural" treatment of unstable particles in 
perturbative field theory is not straightforward. 1 Signal cross sections that are dominated by 
the production and decay of unstable particles with r/m C 1 can be calculated with good 
accuracy in narrow- width approximation (NWA). This and similar approximations, like the 
leading-pole approximation, focus on contributions on or close to resonance and thus greatly 
simplify calculations, since the production and decay of unstable particles (largely) factorises. 
They have been widely employed to predict inclusive and exclusive cross sections. Their use 
to determine background rates for experiments with restrictive selection cuts that eliminate 
resonant contributions and emphasise peripheral phase space regions can be problematic. 
In such cases, users of general-purpose event generators like PYTHIA [[| and HERWIG |3j 
have applied a suggestive procedure that combines results in NWA at the cross section level. 
For ti backgrounds, for example, ti and Wtb results are added to account for double- and 
single-resonant contributions. Since ti calculations in NWA implicitly contain sub- and non- 
resonant contributions that have been integrated out, this procedure can lead to significant 
double-counting |§]. It also neglects interference effects: top pair production and associated 
top (Wtb) production are specialisations of one and the same process, since initial and final 
states of both "processes" are identical. Furthermore, calculations in NWA often do not 
include full spin correlations. 

1 Even limited, appealing schemes like the fermion-loop scheme become rather involved for all but the 
simplest applications || ^| . 
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Examples of new particle searches at present and future hadron colliders with substantial 
tt and Wtb backgrounds include H -> W + W~ 0, g |, [10], 0, and # -> r+r" @, 
13] , [L4| decays, leptonic signals for cascade decays of supersymmetric particles [jn|, the 



Randall- Sundrum radion decays <p ~* W + W~ — > i+i'vv [TJ|] and <p —>■ hh —>■ bbr + r~ |T7 
and searches for if" — > z/ in models with a singlet neutrino in large extra dimensions 
I8| . Reliable phenomenological studies of these searches require tools that allow accurate 



calculations of top pair production and decay in resonant as well as non-resonant phase 
space regions and that are not susceptible to the shortcomings mentioned above. The 
calculations should therefore employ complete matrix elements, i.e. the sum of all leading- 
order (LO) amplitudes. These matrix elements - in fact all resonant fixed-order amplitudes 
- exhibit unphysical singularities, and a finite-width scheme has to be applied to reflect that 
in field theory propagators of unstable particles acquire complex poles when self-energies 
are resummed to all orders. The set of higher-order contributions that has to be included 
to adequately model finite-width effects is not uniquely determined. Moreover, variations of 
higher order in T/m, as well as exclusion of problematic phase space regions, e.g. thresholds, 
are permissible. Consequently, a variety of competing schemes exists |j[ ^ [19], pl|| . 

The purpose of this paper is to compare leading-order tt cross sections calculated in 
narrow-width or double-pole approximation (DPA) to cross sections that take into account 
all sub- and non-resonant amplitude contributions, and to investigate the consistency of sev- 
eral practical finite-width schemes. The program we developed for this purpose is described 
in Section III], with particular emphasis of finite-width schemes and their implementation. 
In Section |III A results for inclusive top pair production are presented, followed by results 
for important top backgrounds to new particle searches in and beyond the Standard Model 
(SM) in Sections 11 IB and III Q , respectively. In Section [V, we conclude with a summary 
and outlook. 



II. PROGRAM DESCRIPTION 

We introduce a LO program for tt production at hadron colliders with up to two additional 
jets that is not specialised to resonant phase space regions and hence has to include complete 
tree-level matrix elements for the contributing 2^6, 2^7 and 2^8 subprocesses. If the 
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W + W decay products are abbreviated as W, tt production includes the subprocesses 

gg^bbW, qq^bbW, (1) 

tt + 1 jet production includes the subprocesses 

gg -> bbW g , qq-^bb\V g , qg ^bbW q, qg ->bbW q, (2) 

and tt + 2 jet production includes the subprocesses 

gg ^ bbW gg , qq->bbW gg , qg^bbW qg , qg ^ bbW qg , 

gg — > bb W gg , qq —> bb W qq , qq ^ bbW qq , qq — > 66 W gg . (3) 

The program contains subprocess matrix elements for the di-lepton, single-lepton and all-jet 
decay modes, or more specifically for the following W + W~ decay final states: 

W di _lcpton = tv Cv , (4) 

w sin gle-lepton — ^ v QdQu ; (5) 

Wan-jet = g«gd gdg« • (6) 

For the di-lepton and all-jet decay modes, the program allows to calculate different-flavor 
samples, e.g. with W = e + z/ e /i - ^, as well as same-flavor samples, e.g. with W = e + z/ e e~u e . 
Additional amplitudes with (7, Z — > x (Z — » ^z/^) and (p, 7, Z — > g u g u ) x (p, 7, Z — > 

g^gd) fragments contribute in the di-lepton and all-jet decay modes, respectively. Moreover, 
a finite- width scheme has to be chosen when formulating complete LO matrix elements with 
unstable particles to avoid unphysical singularities in resonant phase space regions that can 
be removed by including contributions to all orders in perturbation theory. 2 Since no known 
scheme is satisfactory in every respect, a cross section by cross section comparison of several 
schemes with complementary properties is suggestive, but requires more than one version of 
each subprocess matrix element defined above. 

Evidently, the creation of all required matrix elements is a considerable task and calls for 
automation. While the program is generally written in C++ to permit greater code locality 
and expressiveness, we prefer faster Fortran code for the matrix element evaluation, since its 
speed determines the program runtime after initial adaptation. Furthermore, to minimize 



The Dyson resummation of top quark self-energy contributions is described in Ref. [Q . 
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the matrix element code, the program should use helicity amplitudes in unitary gauge that 
neglect CKM mixing. MADGRAPH/HELAS [^, [23| is a matrix element generation sys- 
tem that matches our requirements and has recently been extended to processes with 8-10 
external particles. Its output is used as starting point for the matrix element code in our 
program. 

MADGRAPH/HELAS matrix elements use the fixed-width scheme (FWS). 3 In the FWS, 
all propagators of unstable particles are modified according to the following prescription: 

* ^ 1 . (7) 

p2 _ m 2 p2 _ m 2 _|_ 

This substitution is easy to implement, but the resulting matrix elements with Breit-Wigner 
propagators are not gauge- invariant. As discussed in Refs. M, |5], |5[], calculations that employ 
gauge- variant amplitudes and receive sizable contributions from sensitive phase space regions 
can yield highly erroneous results. To remedy this deficiency, various approaches have 
been suggested in the literature that yield manifestly gauge-invariant matrix elements. The 
theoretically most appealing approach is arguably the fermion-loop scheme [fl ||. We do 
not consider it further here, since it is not applicable to processes with unstable particles 
that decay into bosons, including ti production. Even if it were applicable in the case at 
hand, it would require as prerequisite an analytic calculation of effective vertices that has 
not been automatised yet. Its implementation is therefore not straightforward for complex 
multi-particle processes with several types of unstable particles. For the studies in Section 
PTI| we therefore implement two practical finite-width schemes that allow automatic matrix 
element generation for arbitrary processes and guarantee electroweak and 577(3) gauge- 
invariant results: the complex-mass scheme (CMS) |2D| and the overall-factor scheme (OFS) 
21 . 

The CMS introduces Breit-Wigner propagators in a gauge- invariant manner by replacing 
the masses of all unstable particles with a complex value as follows: 

m — > V m? — imT . (8) 

This substitution is performed unconditionally and yields, for example, for the top propaga- 
tor a different expression than the FWS: + \Jvn? t — im t T t ) / (p 2 — mf + im t T t ) . sin 2 9w and 

3 We use HELAS-3, which implements the fixed-width scheme. Note that the widely-used version 2 of 
HELAS implements step-width Breit-Wigner propagators, i.e. l/(p 2 — m 2 + imT 6(p 2 )). No notable 
deviations occur in general, since \p 2 — m 2 \ 3> mT if p 2 < and T/m <ti 1. 
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dependent quantities also acquire complex values in this scheme, since cos 6w = mw/rriz- 
The CMS matrix elements in our program use HELAS-CMS, a modified version of the 
HELAS library that we created by converting masses and widths from real to complex 
variables. 4 

The OFS conserves gauge-invariance while introducing Breit-Wigner behaviour by mul- 
tiplying the complete LO matrix element (with singular propagators for unstable particles) 
with overall factors: 

2 2 2 2 

p — m p — m 

Jvt-compl. X — 5 5 ■ ; — = Jv\ res BW-prop. ~\~ J^non—res. X — 5 7, ; : 7^ • (9) 

pz _ m z _|_ %m Y pz _ m A _|_ %m Y 

For each unstable particle type, one factor is applied for every time-like momentum combi- 
nation that occurs in propagators of that type. The propagators absorb the corresponding 
factor and transform into Breit-Wigner propagators: 

1 x - . ' . ... (10, 



p2 _ m 2 p2 _ m 2 _|_ ^ m p p2 _ m 2 _|_ 

Amplitudes that are non-resonant with respect to a particular momentum combination do 
not absorb the corresponding factor, as indicated in Eq. @. To facilitate the automatic 
construction of OFS matrix elements a scripting-language program was written that scans 
MADGRAPH output, and analyses the structure of all contributing amplitudes. Potentially 
resonant propagators, where one side is only connected to final state particles, are identified 
and the required overall factors deduced. The script then constructs the overall factor prod- 
uct for each amplitude, and outputs Fortran code that calculates the OFS matrix element. 
To optimise the code, combinations of overall factors that occur multiple times are evaluated 
once and the results are reused. 

The comparisons presented in Section |T| also require the calculation of cross sections 
in double-pole and narrow-width approximation. For that purpose, a second program - in 
nature similar to the one used to generate OFS matrix elements - eliminates all amplitudes 
that do not contain potentially resonant t as well as t propagators, thus extracting all double- 
resonant amplitudes with respect to top decay. The generated Fortran code employs the 



4 Complex widths are introduced since MADGRAPH output uses a real constant ZERO for both, vanishing 
masses and widths, as argument for HELAS calls. In CMS matrix element code we define ZERO as 
complex parameter and then also have to declare all widths as complex variables to be compatible. Note 
that all width variables are set to zero in CMS matrix elements, since the widths are contained in the 
mass variables. 
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fixed-width scheme and is used in our DPA calculations. The DPA matrix elements are also 
used in our NWA calculations. To preserve all spin correlations, we choose to implement 
the NWA directly by calculating with off-shell intermediate top quarks in the T t — > limit. 
To that end, the top width is scaled down to T tcfi = eT t , and \A4\ 2 is multiplied by e 2 to 
restore the proper normalisation of the total amplitude. For one resonant propagator one has 
l-Mefrl 2 = 1/e x \M.\ 2 . A setting of e = 1/1000 is used in the program and yields excellent 
agreement with NWA implementations with on-shell intermediate top states. In DPA or 
NWA mode, the program uses a Breit-Wigner mapping for each resonant top propagator 
that covers a limited range of invariant top quark masses. Neglected contributions from 
outside this range introduce a non-statistical error. For the background calculations in 
Sections |III B| and [III C| , off-shell top masses were generated in m t ± 65T t limiting neglected 



contributions to approximately 1% (see Eq. (19) in Ref. ||]). For the inclusive calculations 
in Section |III A| , we increased the range factor to 6500, reducing this error contribution to 



0.01%. A comparison of results given in Table | below with results in Table II in Ref. Q 
confirms that a range factor of 65 is not sufficient when a total error of less than 1% is 
desired. 

When cross sections for tt production with additional jets are calculated with complete 
matrix elements, one finds that computational complexity increases by a factor of more than 
10 for each additional final state particle beyond the tt level. Resulting program runtimes 
quickly exceed what would be considered acceptable for phenomenological studies. To obtain 
the tijj results presented in Section |111 B[ it was therefore necessary to use state-of-the-art 
integration techniques and to develop a method to distribute the Monte Carlo sampling 
over a larger number of processors. The result is OmniComp, a Monte Carlo integration 



framework based on the adaptive multi-channel techniques introduced in Refs. 0, |25], |26 
that allows to conveniently distribute the calculation over many processors in one or more 
computer clusters. 5 OmniComp further accelerates the computation of hadron collider cross 
sections through adaptive Monte Carlo summation of subprocess <S> helicity combination 
channels. The mapping of sub- and non-resonant phase space regions follows the approach 



laid out in Ref. [27|. OmniComp and the phase space mapping library are described in 



more detail in Ref. [28]. A number of tests were applied to verify the correctness of the 



We successfully ran programs on up to 16 processors. 
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program. First, the Lorentz-invariance of the MADGRAPH-generated FWS matrix elements 
was tested. 6 CMS matrix elements were tested by comparison with corresponding FWS 
matrix elements after the complex masses and widths in HELAS-CMS had been set to 
their usual, real values. The automatic generation of OFS matrix elements was tested by 
comparing with the manually created OFS matrix elements of Ref. H. The DPA/NWA 
matrix elements were verified by comparing NWA cross sections with results from programs 
with on-shell intermediate top quarks. The phase space and PDF integration has been tested 
by comparing with known cross sections for the LHC and Tevatron. Moreover, the addition 
of hadron collider capabilities to the general purpose packages O'Mega & Whizard p^] and 
AMEGIC++ H reached the final stage this year, and a comparison of top production cross 
sections to cross-check our implementations is planned for the near future. 



III. NUMERICAL RESULTS 



In this section, we use the program described above to study the difference between cross 
sections with on-shell (NWA) and off-shell (DPA) intermediate top quarks, to determine 
the size of corrections when complete LO matrix elements are included, and to search for 
deviations between results obtained with different finite- width schemes. 7 We first investigate 
these issues for inclusive top pair production and then turn to suppressed top backgrounds 
that are important for new particle searches in and beyond the Standard Model. 

To cover the energy range of existing and future hadron colliders, cross sections are 
calculated for the Tevatron (pp, y/s = 2 GeV), the LHC (pp, y/s = 14 TeV) and a Stage- 
1 VLHC (pp, y/s = 40 TeV). Unless otherwise noted, all calculations use the following 
parameters: m z = 91.187 GeV, G F = 1.16637 x 1(T 5 GeV" 2 and a(m z ) = 1/128.92, which 
translates at tree-level to sin 2 9w = 0.23105 and mw = 79.9617 GeV, as well as the masses 
m t = 175 GeV, m& = 4.4 GeV and m# = 115 GeV. LO formulas for the decay widths then 
yield T t = 1.56 GeV, = 2.01 GeV, Y z = 2.42 GeV and T H = 0.00323 GeV. CTEQ6L 
LO parton distribution functions are employed by default, with a s (mz) = 0.118 and the 
NLO formula. Factorisation and renormalisation scales are fixed at the top mass, except 

6 The Lorentz-invariance of one and the equivalence of two matrix element routines was tested as described 
in Ref. (§. 

7 Note that all results calculated with our program include full spin correlations (see Section ||) . 
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for studies where cross sections with additional jets are taken into account (e.g. in Table 



II|). In this case, the factorisation scale is chosen as \Xf = min(mT) of the top quarks and 
additional jets. This factorisation scale definition avoids double-counting of contributions 
that have already been integrated out in the PDFs. The overall strong coupling constant 
factor is calculated as (a s ) n = n™ =1 a s (m T j), again using the transverse masses of both top 
quarks and any additional jets as input. ATLAS detector resolution and b decay effects 



are modelled as described in Ref. but tagging efficiencies are not taken into account. 
Monte Carlo integration errors are 0.1% or less for inclusive cross sections and 1% or less 
for background cross sections. 



A. Inclusive Production 

For Table |, we choose inclusive ti production and decay into the di-lepton final state 
bb e + v e /i - ^, and compare cross sections in NWA to cross sections with complete matrix 
elements. First, the size of changes is compared for LHC and Tevatron collisions. Generally 
one would expect finite-width effects to be of order T t /m t = 0.009, and the Tevatron correc- 
tion is indeed less than 1%. The LHC cross sections, however, is enhanced by a significantly 
larger factor of 1.06. This effect is not caused by averaged v. exact spin correlations, since 
our NWA results include full spin correlations. Furthermore, less than 1% of the increase 
can be attributed to double-resonant off-shell effects (as seen in Table JJ). The increase is 
mainly caused by previously omitted sub-resonant contributions. These contributions are 
also included for the Tevatron. That no sizable increase occurs there can be traced to the 
fact that the Tevatron cross section is dominated by quark scattering, while the LHC cross 
section is dominated by gluon scattering. As shown in the two rightmost columns, the large 
increase is specific to the gluon- initiated process. In fact, the qq cross section is slightly 
reduced when sub- and non-resonant amplitudes are included. A suggestive kinematical 
interpretation that relates the large increase to off-shell contributions at hard scattering en- 
ergies below the on-shell top pair production threshold that are amplified by steeply falling 
PDFs is therefore misleading. A comparison of calculations with CTEQ4L and CTEQ6L 
shows that recent PDF improvements decrease cross sections uniformly by 18%, thus having 
little effect on the ct fws /cxnwa enhancement factor. For inclusive cross sections, differences 
between finite- width schemes are expected to be of higher order in T t /m t . The third com- 
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TABLE I: Cross sections in NWA and with complete matrix elements for inclusive ti production 
and di-lepton decay (bb e + u e Effects are compared for colliders, PDF sets and finite-width 

schemes. Tevatron corrections are of order Tt/mt = 0.009, but LHC corrections are larger. PDF 
improvements decreased the LHC cross section by 18%, but the correction is robust. FWS, CMS 
and OFS yield consistent results. All cross sections are given in fb. Note that NWA cross sections 
include full spin correlations (see Section [n|). 



Collider 


CNWA 


CFWS 


°"FWS 
"NWA 


CT 99,FWS /f 99 ,FWS\ 
""SS.NWA ^ "FWS ' 


°"ijf,FWS /Oqq.FWSN 
Cqq.NWA V ""FWS ' 


LHC 


5.86 x 10 3 


6.19 x 10 3 


1.06 


1.06 (88%) 


1.00 (12%) 


Tevatron 


66.5 


66.5 


1.00 


1.10 (5%) 


0.99 (95%) 



PDF set 


cnwa 


0"FWS 


<7FWS/C"NWA 


CTEQ6L 
CTEQ4L 


5.86 x 10 3 
7.18 x 10 3 


6.19 x 10 3 
7.58 x 10 3 


1.06 
1.06 



Scheme 


a 


CX/<7NWA 


NWA 


5.86 x 10 3 


1.00 


FWS 


6.19 x 10 3 


1.06 


CMS 


6.19 x 10 3 


1.06 


OFS 


6.18 x 10 3 


1.06 



parison in Table | shows that FWS, CMS and OFS yield results that agree when integration 
errors of 0.1% are taken into account. 

Cross section changes when progressing from NWA to DPA and finally to complete LO 
matrix elements are displayed in Table JT] for the di-lepton, single-lepton and all-jet channels 
of inclusive ti production at three hadron colliders: Tevatron, LHC, and a Stage-1 VLHC 
with yfs = 40 TeV. In all cases the ratios (Jcms/^nwa and ctcms/cdpa are very similar: 
Replacing on-shell with off-shell intermediate top quarks changes cross sections by no more 
than 1%. Effects are generally small at the Tevatron. At pp supercolliders, on the other 
hand, cross sections increase uniformly by about 7% for all decay modes and collider energies. 
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TABLE II: Change of cross sections with complete matrix elements relative to cross sections in 
NWA for inclusive ti production and di-lepton, single-lepton and all-jet decay at the Tevatron, 
LHC, and Stage- 1 VLHC (y/s = 40 TeV). At pp supercolliders, cross sections increase uniformly 
by about 7% for all decay modes and collider energies. Effects are similar for cross sections in DPA 
(shown in parenthesis): replacing on-shell with off-shell intermediate top quarks changes cross 
sections by 1% or less. The cuts pt > 15 GeV, \r]\ < 4.5 and AR > 0.6 are applied for channels 
with singular phase space regions related to massless particles. 





ccms/cnwa (ocms/cdpa) 


di-lepton 


single-lepton 


all- jet 


Tevatron 


1.00 (1.01) 


1.01 (1.02) 


1.00 (1.01) 


LHC 


1.06 (1.06) 


1.07 (1.07) 


1.07 (1.07) 


VLHC 


1.06 (1.06) 


1.07 (1.06) 


1.08 (1.07) 



B. Backgrounds to SM Higgs Searches 

The unexpectedly large increase of inclusive cross sections at pp supercolliders when sub- 
and non-resonant contributions are included raises the question how large corresponding ef- 
fects are for new particle searches with significant ti backgrounds. When optimised selection 
cuts are used to suppress top pair production, the dominant double-resonant it contributions 
are typically suppressed by factors of order 10~ 4 and the importance of contributions from 
sub- and non-resonant phase space regions can increase considerably. A central jet veto, e.g. 



PTj > 15 GeV and \r]j\ < 3.2 , (11) 

is very effective in suppressing the ti background to inclusive H — > WW searches at the 
LHC H, H, H9- The ve t° °f Eq. ([□]) reduces the inclusive cross section of 6 pb to 14 
fb when the NWA or DPA is applied, whereas a calculation with complete matrix elements 
yields 26 fb. Sub- and non-resonant contributions increase the result by a factor 1.8. As 
in the inclusive case, moving from NWA to DPA or switching finite-width schemes changes 
the corresponding result very little in comparison. The distributions in Fig. [I] show this 
relationship for differential cross sections. 
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FIG. 1: Charged lepton separation in rj-cp space and transverse mass distributions for suppressed tt 
production and di-lepton decay (bb e + v e H'v^) at the LHC. The applied central jet veto (pr > 15 
GeV and \rj\ < 3.2) reduces the acceptance to 4 x 10~ 3 . Differential cross sections with complete 
matrix elements (solid lines) and in NWA and DPA (dashed lines) are shown. Sub- and non- 
resonant amplitude contributions enhance the total cross section by a factor 1.8. Off-shell top effects 
(DPA v. NWA) and deviations between finite-width schemes (FWS, CMS, OFS) are negligible. The 



transverse mass is defined as = \/2p^$J T [l — cos A9(££, $ T )]. 



TABLE III: Top background cross sections with up to two additional jets for the H — > WW — > 
e^n^ftj, decay search in weak boson fusion at the LHC. The light Higgs-optimised selection cuts 
and event classification from Ref. |llj are applied. All cross sections are given in fb. Sub- and 
non-resonant amplitude contributions enhance the total top background by a factor of 1.1. The 
tt background without additional jets is strongly suppressed because in this case both b quarks 
need to be resolved as forward jets with wide separation in pseudorapidity and very large di-jet 
invariant mass. 





tt 


ttj 


tijj 




a 


c/cNWA 


a 


c/o"NWA 


a 




NWA 


0.020 


1.0 


0.94 


1.0 


0.24 


1.0 


FWS 


0.044 


2.1 


1.08 


1.1 


0.24 


1.0 


CMS 


0.044 


2.1 


1.07 


1.1 


0.24 


1.0 


OFS 


0.044 


2.1 


1.07 


1.1 


0.24 


1.0 
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H — > WW searches are usually tuned for intermediate Higgs masses around 170 GeV, 



where the H — > WW branching ratio is large. In Ref . |TTJ , the search for H — > WW decays 
in weak boson fusion at the LHC is studied for the light Higgs (m# = 115 GeV) favoured 
by LEP experiments JHJ]. The additional forward jets in weak boson fusion permit powerful 
selection cut optimisations that make this search channel competitive - a 5a discovery is 
possible with 35 fb _1 - even for relatively low Higgs masses where the H — > WW branching 
ratio is small. In this search scenario, ti + jets production is the dominant background and 
its accurate determination is essential. The ti background is strongly suppressed, because 
for final states without additional jets both b quarks need to be resolved as forward jets with 
wide separation in pseudorapidity and very large di-jet invariant mass. In Ref. [11]], complete 
matrix element corrections were calculated for the ti and tij backgrounds using the OFS. As 



shown in Section [HI Q , OFS cross sections can be artificially reduced. Our program allows 
to calculate these corrections using the FWS, CMS or OFS. The results are given in Table 



1 I J and show that the OFS is reliable in this case. Table [HI] also displays first results for 
ti+2 jets production calculated with complete LO matrix elements in the literature. Sub- 
and non-resonant amplitude contributions enhance the total top background by a factor of 
1.1. 



C. Backgrounds to Beyond-SM Physics Searches 

At the LHC, top backgrounds also play an important part in searches for physics beyond 
the Standard Model. In this section we focus on two studies where ti production constitutes 
the dominant background: the search for hadronic r decay of a heavy charged Higgs in 
supersymmetric models and the radion decay <fi — > hh — > bbr + r~, where one r decays 
leptonically and the other hadronically. 

The production of a charged Higgs with m#± > m t in supersymmetric models at high 
tan/? was analysed in Ref. p2] . Production proceeds through gg — ► H^tb and is followed 
by the decays H ± — > rv (with hadronic r decay) and t — > jjb. Applying the selection cuts 
of this ATLAS study, we calculate ti background cross sections in NWA and with complete 
matrix elements to determine the enhancement factor. The results are shown in Table [TV]. 
Sub- and non-resonant amplitude contributions enhance the top background by a factor of 
1.1. The ATLAS analysis takes sub- resonant contributions into account by combining ti and 
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TABLE IV: Top background for heavy charged Higgs production gg — > H^tb and decays — > rz/ 
(with hadronic r decay) and i — > jjb at the LHC. The selection cuts of the ATLAS analysis in 
Ref. are applied. Sub- and non-resonant amplitude contributions enhance the top background 
by a factor of 1.1. The ATLAS analysis combines ti and Wtb results in NWA at the cross- 
section level. This procedure can lead to substantial double-counting of sub- and non-resonant 
contributions 0, evidently a 30% effect in the case at hand. All cross sections are given in fb. 
Parton-level results are rescaled by a factor 0.16, so that our NWA result and the PYTHIA- 



ATLFAST P, B3 tt result given in Table 3 in Ref. [32] match. 





a 


d/(T NW A 


NWA 


0.343 


1.00 


NWA (ti + Wtb) a 


0.485 


1.41 


FWS 


0.376 


1.09 


CMS 


0.378 


1.10 


OFS 


0.364 


1.06 



"calculated in Ref. 32 



Wtb results in NWA at the cross-section level. This procedure can lead to substantial double- 
counting of sub- and non- resonant contributions [§J . Our enhancement factor indicates that 



the actual top background is 23% lower than the estimate in Ref. f32j. In Ref. the 
analysis was extended to the search for H — > t^v in models with a singlet neutrino in 
large extra dimensions, and we expect a similarly reduced top background if sub-resonant 
contributions are included at the amplitude-level. 

In Table [V|, enhancement factors are given for the top background to the decay — > 
hh — > bbr + T~ of a Randall-Sundrum radion with mass 300 GeV. The two r leptons decay 
leptonically and hadronically, respectively. The selection cuts of the ATLAS analysis in 
Ref. |T7| are applied. Specific model parameters are given in the table caption. The effect 
of sub- and non-resonant amplitude contributions is small for this top background, in fact 
smaller than for inclusive ti production. 

The OFS cross sections in Tables [IV] and are several percent lower than the corre- 
sponding FWS and CMS cross sections, which agree within the integration error of 1%. To 
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TABLE V: Top background for the radion decay <f> — ► hh — ► bbr + T~ at the LHC, where one r decays 
leptonically and the other hadronically. The selection cuts of the ATLAS analysis in Ref. |l7j are 
applied. The radion vacuum expectation value = 1 TeV, the radion-SM Higgs mixing parameter 
£ = 0, the radion mass = 300 GeV and the lightest Higgs mass = 125 GeV. Our results 
indicate that sub- and non-resonant amplitude contributions change the top background by not 
more than 3%. All cross sections are given in fb. Parton- level results are rescaled, so that our 
NWA result and the PYTHIA-ATLFAST g ||] it result given in Table 5 in Ref. @ match. 





a 


c/onwa 


NWA 


3.27 


1.00 


FWS 


3.36 


1.03 


CMS 


3.34 


1.02 


OFS 


3.17 


0.97 



understand why the OFS may not be suitable for all ti background calculations, we integrate 
the small phase space region where both top quarks and the intermediate Z boson are close 
to resonance. More specifically, we require \mwb — mt\ < 2r 4 and \mww — m z\ < ^z- The 
results are shown in Table [V|. The DPA is excellent and contributions from amplitudes with 
resonant Z propagator are negligible. In the OFS, the dominant double- resonant top ampli- 
tudes are artificially suppressed by the overall factor — m z )/(p z — m z + im z T z )\ 1. 
The resulting OFS cross section is consequently much smaller than CMS or FWS cross 
sections. This example illustrates that cross sections for multi-resonant processes cannot 
be calculated reliably with the OFS if sizable contributions arise from phase space regions 
where several amplitudes with different resonance structure compete. Artificially reduced 
cross sections can even occur for single-resonant processes, given that non-resonant contri- 
butions are sizable in phase space regions close to resonance. The authors of Ref. B, for 
example, compared OFS and fermion-loop scheme cross sections for radiative W production 
and found that OFS results are 30% lower close to threshold. Despite some effort 8 , we were 
unable to find a phase space region, where CMS and FWS cross sections for ti production 

8 We considered, for example, the selections m\yb > 500 GeV, \rfb\ > 3.5 and ARww < 1-0 for collider 
energies up to 100 PeV. 
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TABLE VI: Cross sections for tt production with di-lepton decay (bb e + v e n~Va) at the LHC, 
calculated in DPA and with complete matrix elements using several finite-width schemes. Only 
the phase space region where \ mww— m z\ < and \ mwb— m t\ < 2I\ is integrated. In this region 
Z boson and top quark propagators are resonant. The DPA is excellent, i.e. resonant tt production 
dominates. The contribution from amplitudes with a resonant Z propagator is negligible. In this 
phase space region OFS matrix elements are inadequate, since the dominant double-resonant top 
amplitudes are artificially suppressed by the overall factor \(p 2 z — m 2 z )/(p z — m 2 z + imzTz)\ "C 1. 
The OFS result is 30% smaller than the CMS and FWS results. All results are given in fb. 





a 


0"/o"DPA 


DPA 


0.0168 


1.00 


FWS 


0.0170 


1.01 


CMS 


0.0170 


1.01 


OFS 


0.0118 


0.70 



showed significant discrepancies. We therefore conjecture that calculations employing the 
gauge-variant fixed-width scheme may be used to obtain reliable predictions for the pro- 
cesses considered here. We note that the reliability of the fixed-width scheme has recently 
also been established for e + e~ — > 6 fermion processes |H . 



Our LO calculations do not include log-enhanced higher-order contributions from collinear 
g — > bb configurations for initial state gluons. An improved treatment would include gb 
scattering matrix elements convoluted with the b quark PDF. Then, a subtraction of the 
gluon splitting term would also be required to avoid double-counting f35| . However, the 
additional net contribution to inclusive tt production is less than 2% [[JIJ and can safely be 
neglected in our analysis in Section [111 A| . In the weak boson fusion Higgs search discussed 
in Section flIB|, one or both b quarks have no finite transverse momentum threshold, but 
collinear contributions are small within typical cuts j§]. For central jet veto suppressed top 
backgrounds, on the other hand, one would expect more pronounced collinear enhancement. 
The selection cuts for the tt background studies in Section |III Q require that the b quarks 
are resolved with a transverse momentum of at least 15 GeV. The collinear region is thus 
avoided. 
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IV. SUMMARY AND OUTLOOK 



We presented cross sections for top quark pair production at hadron colliders with up 
to two additional jets resulting in 6-, 7- and 8-particle final states calculated with complete 
tree-level matrix elements. Our program includes di-lepton, single-lepton and all-jet decay 
modes, and implements three practical finite-width schemes, i.e. the fixed-width, complex- 
mass and overall-factor schemes, as well as the narrow-width and double-pole approximation 
for comparison. While our LO calculations are subject to substantial scale uncertainties, the 
obtained cross section ratios are expected to be robust. For inclusive production, advancing 
from NWA to DPA by replacing on-shell with off-shell intermediate top quarks, changes cross 
sections by 1% or less. The inclusion of sub- and non- resonant amplitudes increases NWA 
or DPA cross sections by 5-8% at the LHC and VLHC, but has little effect on Tevatron 
cross sections. Top backgrounds to new particle searches are often suppressed by optimised 
selection cuts that can enhance the importance of sub- and non-resonant contributions con- 
siderably. We updated on-shell ti/Wtb background estimates for the H — > WW decay in 
weak boson fusion, the hadronic r decay of a heavy and the <fi —>■ hh — > rrbb radion 
decay at the LHC, and found corrections from 3% to 30%. All calculated FWS and CMS 
cross sections agree within errors. Gauge- violating effects of the FWS appear to be generally 
negligible for the processes considered here. Our calculations show further that the OFS 
may yield underestimated cross sections and should be applied with caution in studies with 
suppressed top backgrounds. 

Because of the large scale uncertainties of LO cross sections, precise absolute predictions 
for top pair production at hadron colliders cannot be achieved with tree-level calculations. 
The extension of LO to NLO calculations in the framework of the narrow-width and double- 
pole approximations was first explored in the context of weak gauge boson production |3"7[ 
and has recently also been carried out for top pair production at hadron colliders ||38f| . The 
results in Section |III A| imply that sub-resonant contributions need to be included in NLO 
calculations for inclusive ti production at pp supercolliders to achieve a theoretical error 
of 0(5%). A common method to improve LO predictions for suppressed top backgrounds 
is to apply a reaction-specific .fT-factor, i.e. to rescale all LO results by multiplying with 
K = <Ji nc i, nlo/ >cr ind, LO- When sub-resonant and non-resonant phase space regions contribute 
substantially to cross sections, the merit of such procedures has to be tested by comparing 
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with fully differential NLO calculations that cover resonant and non-resonant phase space 
regions. The starting point for a complete NLO calculation of top pair production, i.e. 
a calculation that is not specialised to the double-resonant phase space region, would be 
the evaluation of the NLO corrections of the complete matrix element for the bbW + W~ 
final state. The calculation of the real emission corrections is straightforward, since the W 
bosons are on-shell. However, the evaluation of the virtual corrections of this 2 — > 4 process 
involves 1-loop hexagon amplitudes, whose computation is still very challenging . While 
the tt+ jets program described in this paper allows to calculate the real emission component 
of a complete calculation of pp — > 6 fermions at NLO in QCD, the evaluation of the virtual 
corrections for such 2 — > 6 processes is well beyond present capabilities. NLO predictions for 
many-particle processes with multiple scales can be further improved by resumming higher- 
order contributions with large logarithms, such as a s log(m^/p^ -) in the case at hand. 

In addition to precise and accurate calculations for hard scattering subprocesses, a reli- 
able comparison of theoretical predictions and experimental data also requires the proper 
inclusion of parton showering, hadronisation and detector effects. To standardise the co- 
operation of parton-level Monte Carlo programs (with full matrix elements) and showering 
and hadronisation event generators - which in turn produce input for detector simulations 
- a generic interface has been specified recently in Ref. |4"D |. In the near future, we plan 



to implement this interface and to make our complete LO top pair production program 
available to interested experimental physicists. 
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